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Abstract 

The event by event fluctuations of the ratio of positively to negatively charged 
hadrons are predicted within the UrQMD model. Corrections for finite ac- 
ceptance and finite net charge are derived. These corrections are relevant to 
compare experimental data and transport model results to previous predic- 
tions. The calculated fluctuations at RHIC and SPS energies are shown to 
be compatible with a hadron gas. Thus, deviating by a factor of 3 from the 
predictions for a thermalized quark-gluon plasma. 



^ Feodor Lynen Fellow of the Alexander v. Humboldt Foundation 
E-mail: bleicher@nta21bl.gov, SJeon@lbl.gov, vkoch@lbl.gov 

1 



In this note we exploit a novel idea of a signal proposed by [|IJ based on event by event fluc- 
tuations of the N + /N_ ratio (iV + and iV_ being positively and negatively charged hadrons). 
In particular, we would like to precisely specify the observable which is to be compared with 
the prediction of Ref. 0. The magnitude of these fluctuations was estimated from thermo- 
dynamics to be a factor of 3 different for a QGP compared to a conventional hadronic gas. 
The predictions of ]TJ were derived using a grand canonical ensemble and assuming a zero 
net charge. In a grand canonical ensemble charge is conserved only on the average, whereas 
charge is conserved event-by-event in real heavy ion collisions and the model calculations 
presented here. The assumption of a grand canonical ensemble thus is only valid in the limit 
of P = (N ch ) Ay / (N ch ) total — > and correction to that are of the type (1 — P) |j]|§- Thus, 
in order to compare experimental data with the predictions of [|TJ, corrections have to be 
applied, which will be presented here. 

For our study of this observable, we employ the Ultra-relativistic Quantum Molecular 
Dynamics model (UrQMD) ||. This model has been shown to reproduce the measured data 
on event-by-event fluctuations in Pb+Pb collisions at the SPS Q. UrQMD is a microscopic 
transport approach based on the covariant propagation of constituent quarks and diquarks 
accompanied by mesonic and baryonic degrees of freedom. It simulates multiple interactions 
of ingoing and newly produced particles, the excitation and fragmentation of color strings 
and the formation and decay of hadronic resonances. At SPS and RHIC energies, the 
treatment of subhadronic degrees of freedom is of major importance. In the UrQMD model, 
these degrees of freedom enter via the introduction of a formation time for hadrons produced 
in the fragmentation of strings |J. The leading hadrons of the fragmenting strings contain 
the valence-quarks of the original excited hadron. In UrQMD they are allowed to interact 
even during their formation time, with a reduced cross section defined by the additive quark 
model, thus accounting for the original valence quarks contained in that hadron ||. The 
thermodynamics properties of the model in equilibrium are similar to a Hagedorn gas with 
a limiting temperature of 135 MeV |J. 

Let us shortly review the main ideas of the charge fluctuations as elaborated in |1] . The 
following connection between the observable charge ratio fluctuations and the fluctuations 
of the elementary net charges of the system has been proposed: 
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D=(N ch )(5R 2 )^A(5Q 2 )/(N ch ) , (1) 

with Q = J2aQaN a , where a runs over all hadron/parton species in the interval Ay <C 2/totai 
in an event and q a is the charge of the species a and N a is the number of particles of species 
a. The ratio fluctuations (SR 2 ) are defined by 

(5R 2 ) = (R 2 ) - (R) 2 , (2) 

where the (...) denotes event averaging. The ratio R = N + /N~ is obtained for each event 
separately in the given rapidity interval Ay by counting the charges of the hadrons. Note 
that weak decays are not taken into account - in line with the experimental setups. 

As already mentioned at the beginning of this letter, the predictions for D obtained in |T]| 
are based on a grand canonical ensemble assuming zero net charge in the system. In reality, 
the system under consideration has a small but finite net charge. Also, the assumption of a 
grand canonical ensemble is only justified in the limit of P —>■ 0. 

For a finite acceptance, transport model calculations and the experimental data need 
to be corrected in order to be compared with the predictions of Jl| (see Eq. 11]). One 
observation is that as P — > 1, the charge fluctuation has to vanish because the global charge 
does not fluctuate. Also Eq.(|I]) is no longer valid in the same limit if the net charge of the 
system under consideration is non-zero. 

The complete analysis is fully explained in the upcoming paper M. In this letter, we list 
the key results. In terms of iV + and AL, the ratio fluctuation is given by || 

where (. . .) A is the event average in a given rapidity window Ay and 

RAy = (N+) Ay /(N„) Ay . (4) 
In an ideal detector with full coverage, the ratio fluctuation is 

( 5jR2 )total = W H ^total TTTT TTTT2 ( 5 ) 

\ iV -/total \ iV +/total 

where Z is the fixed total charge in the system, wh = (5N*) / (NJ) ~ 1, Here (. . -) total 
denotes the event average in the full phase space. At SPS, (<5-R 2 ) total ~ 1% and at RHIC, 



(<5-R 2 ) total ~ 1CT 5 . This is a consequence of the charge conservation. No detailed information 
on dynamics can be obtained from Eq.(|5[). 

Having a finite rapidity window modifies the above full phase space result in an essen- 
tial way. The particle number in a given rapidity bin now has two independent source of 
fluctuations: 

(i) Event-by-event fluctuation in the overall number N r of a given particle species r due 
to different impact parameters and inelasticities. 

(ii) Given a class of events with the same N r , additional fluctuations in a specific rapidity 
window occur due to different distributions of the N r particles in phase space. 

These two fluctuations are independent of each other. The squares of the fluctuations are 
additive. Note that the grand canonical ensemble approach requires Ay <C ytotab therefore 
the fluctuations from (i) are negligible as shown below. 

Assuming that the fluctuations (ii) follow a binomial distribution and that the probability 
P is the same for all species of particles, the number fluctuations in a given rapidity window 
are 

(5N*) Ay = P*(5N*) total + P(l-P)(N) total 

= P 2 m 2 ) totiil + (l-P)(N) Ay , (6) 

where (. . .) A denotes the event average within Ay and we used (N c ^) Ay = -P(-/Vch) t otar 

In the ratio (5N 2 ) A / (N±) 2 Ay needed to construct (5R 2 ) Ay (c.f. Eq.(D), the P 2 factor in 
the first term cancels. (The same holds for the cross terms.) Hence, the contribution of the 
first term is (SR 2 ) totSkl in Eq.fl5|). The second term replaces (R) factor in Eq.(lO) in Ref. || 
to yield 

(5R 2 ) Ay = (SR\ ota , [ JM + (1 - P)R\ y D 2 + _ , (7) 

V^Hotal / 

where -D+_ is the expression given in Eq.(13) in Ref. ]§[]. 

The prediction given in Ref. [Q is for (N^) Ay D\_. To extract D 2 __ from Eq.(|7]), we need 
to apply two correction factors, = R Ay and C y = (1 — P). Thus, collecting the correction 
factors yields 



^-^ - {^ch) A yL>+- + ^ _ p j ^ V ch) tota i— ^ J ~ \^ch) Ay lJ + ^ , (OJ 

where {N c h) A = -P(-/Vch) t otai * s used. For P < 1, the second term is negligible at SPS and 
at RHIC. 

Let us summarize the correction factors that need to be applied to take care of the effects 
of the finite net charge and the finite acceptance window. 

1. In order to correct for the finite net charge within the acceptance due to baryon 
stopping, one has to apply a factor given by 

^ = Rly = 0) 



(JV_> 



Ay 



to the experimental data and the model calculations to compare with the pion gas and 
quark gas result of [[!]. 

2. In order to correct for the finite bin size in rapidity, and in order to incorporate global 
charge conservation one has to rescale the experimental data and the transport model 
predictions by a factor of 

C y = l-P = l- . (10) 

\ iv ch/total 

Thus, the basic observable to be compared with the predictions calculated in |L| is 

(N ch ) Ay (5R 2 ) Ay 



D 



Cy 



1 quark gluon gas 

2.8 resonance gas (11) 
4 uncorrelated pion gas 



Again the subscript Ay denotes the average taken in the rapidity acceptance, while the 
subscript 'total' indicates the average of 47r acceptance. 

Fig. |l] shows the D values as predicted by the UrQMD model (symbols) as compared 
to the estimates for the resonance gas (dashed line) and the quark gluon gas (full line). 
The shown acceptance cuts are chosen according to the experimentally accessible rapidity 
windows of the NA49 (2.5 < y < 4.5) and STAR (—1 < y < 1) detectors. One observes, 
that the UrQMD model predictions of the charge ratio fluctuations are in agreement with 



the expectations of a hadron gas. The predicted fluctuations are a factor of 3 larger than 
the fluctuations expected if a QGP was formed. Thus, the proposed observable - the charge 
ratio fluctuations - predicted for a QGP can not be mimicked by either multiple string 
break-up nor strong rescattering effects which may drive the system towards an equilibrated 
Hagedorn gas behavior. 

Fig. |2] explores the fluctuation parameter D as a function of the width of the inspected 
rapidity window (y cm ± -^) in Au+Au, b < 2 fm at = 200 AGeV. Full squares denote 
the charge ratio fluctuation values obtained with all corrections included as discussed above. 
Open squares show the charge ratio fluctuation values without the correction for finite 
rapidity window and net charge. In line with the findings of ||, one observes a strong 
decrease in the fluctuation values if no correction is applied.^] However, the inclusion of the 
necessary corrections (see Eqs. § and [[(]) yields fluctuations similar to the ones obtained 
from a resonance gas over all inspected rapidities. Notice further the increase of the values 
obtained from UrQMD in rapidity windows of Ay < 2. This increase is due to the vanishing 
correlations from resonance decays which is present at larger rapidity windows. Details are 
discussed in || (see also 0). 

In conclusion, the charge ratio fluctuations at SPS and RHIC energies are predicted 
from the UrQMD model and compared to a hadron gas and QGP estimate. It is shown 
that the UrQMD model results are compatible with the formation of a hadron gas at SPS 
and RHIC energies. The transport model simulations predict fluctuations that are by a 
factor of 3 larger then the fluctuations characteristic for QGP formation. The dependence 
of the fluctuation on the rapidity width is predicted and shown to be compatible with a 
resonance gas if all corrections are included. This observable can be easily studied as 'Year 
1' observable with STAR at RHIC. It can also be directly accessed by the NA49 experiment 
at the CERN-SPS. 

Finally let us stress that it is the corrected observable D defined in Eq. flTT| ) that must 
be compared with our predictions. A measurement of D ~ 1 indicates the presence of a 



1 In Ref. |J the authors used PYTHIA to calculate (Y ch ) Ay (5R 2 ) Ay . Here we are using UrQMD. 
Nonetheless, we get identical results as shown in Fig. |2| 
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QGP state in the system created by the heavy ion collision. 
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FIG. 1. The corrected fluctuation parameter D = (N^) Ay(8R 2 ) Ay/(Cu,Cy) as a function of 
centrality. Circles denote Pb+Pb collisions at 160AGeV, while squares denote Au+Au collision 
at yfs = 200 AGeV as predicted by the UrQMD model. The dashed line show the analytical 
predictions for a hadron gas pj, while the full line depicts the value obtained from lattice QCD 
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FIG. 2. Fluctuation parameter D = (N ch ) Ay {5R 2 ) A y/(C^Cy) and D = (N ch ) Ay (5R 2 ) \y as a 
function of the width of the rapidity window (y cm ± ^) in Au+Au, b < 2 fm at = 200 AGeV 
as predicted by the UrQMD model. Full squares denote the charge fluctuation values D obtained 
with all corrections included. Open squares show the charge fluctuation values D without the 
correction for finite rapidity window and net charge. The lines are the same as in Fig. [l]. 
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